Periodic orbit quantization of a Hamiltonian map on the sphere 
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I. INTRODUCTION 

The general theory of periodic orbit quantization was first developed many years ago by 
Gutzwiller [|1], |^ and has ever since played an important role in our understanding of chaotic 
systems in the semiclassical regime. By applying stationary phase approximations to the 
trace of the Green's operator, Gutzwiller was able to derive a semiclassical formula for the 
density of states as the sum of a smooth part and an oscillating sum over all periodic orbits 
of the corresponding classical system. Gutzwiller 's 'trace formula', however, is only valid for 
isolated periodic orbits which, in general, restricts his results to fully chaotic systems where 
the dynamics is completely hyperbolic. 

At the other end of the scale lie integrable systems. Armed with the knowledge of all 
constants of motion one may obtain the semiclassical eigenenergies using EBK torus quanti- 
zation 0. Alternatively, Berry and Tabor [^] showed that EBK quantization could be recast 
as a 'topological sum' over the periodic orbits and obtained an analogue of Gutzwiller's trace 
formula for integrable systems. The Berry- Tabor formula has also been extended into the 
near-integrable regime p, |^. However most physical systems are of a mixed type, exhibiting 
both regular and chaotic behaviour. Some progress on the periodic orbit quantization of 
mixed systems has also been made [0, H, 0. 

The theory referred to above is for autonomous systems. The dynamics of periodically 
driven systems and quantum maps is best described by the so-called Floquet operator [jlO| . 



Upon deriving an approximation for the trace of this unitary operator in terms of the classical 
periodic orbits, one may obtain a semiclassical approximation for the spectral density of its 
eigenphases. This trace formula has been derived for a class of maps on the plane ||1T|, Il2|, 



on the torus |]T3| and for an array of special mappings with mixed results. The trace formula 
for Arnold's cat map is exact |T^, [1^. However the more general sawtooth map gives poor 
results [|13|, |T6|, It was found necessary to include boundary contributions in the trace 
formula of the baker's map ||18l 119], EH, |2l|, |22[, while ghost contributions needed to be 



included in the case of the kicked top ^ ^ and the standard map ^ ^ 

Poincare maps have also been considered 



The model for which we will apply Gutzwiller's method of periodic orbit quantization 
is quite dissimilar to the above mappings. The dynamics of our model is dominated by an 
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infinite number of isolated stable periodic orbits of arbitrarily long period. The resonances 
associated with these periodic orbits are circular disks and fill in the phase space in a manner 
resembling circle packings. Three of these mappings were introduced in a previous paper |^ 



each corresponding to a particular phase space geometry: planar, hyperbolic and spherical. 
We will consider only the spherical map in this paper. Its quantum analogue lives in a finite 
dimensional Hilbert space allowing us to circumvent convergence problems associated with 
semiclassical trace formulae |^ . In addition, our model does not suffer from the exponential 
proliferation of periodic orbits found in most chaotic systems. The task of finding all periodic 
orbits of a given period is much simplified with the help of symbolic labeling and the use of 
special symmetries. Our mapping is not chaotic in a strict sense. However the presence of 
infinitely many stable periodic orbits with arbitrarily long period generates highly irregular 
motion. This behaviour is located on the borderline between regular and chaotic motion, 
but in a manner very different from the generic case of KAM systems. The strong deviation 
is explained by the mapping's non-smooth construction. 

Our motivation for considering such a system is as follows. Tests of the applicability of 
periodic orbit quantization as a theory which adequately describes quantum systems in their 
semiclassical limit has, for the most part, been focussed on fully chaotic systems. Systems 
with regions of stability have been ignored due to complications in the theory which arise 
whenever non- isolated periodic orbits are present. For integrable systems, the presence of 
periodic tori surrounding each stable periodic orbit is accommodated in the theory of Berry 
and Tabor. However when such a system is perturbed into the near-integrable regime, 
these tori resonate and form island chains of high-order stable and unstable periodic orbits 
of equal period. If the perturbation is small, the newly created periodic orbits will not be 
sufficiently isolated, and hence, Gutzwiller's theory cannot be applied. In such cases, uniform 
approximations need to be made which bridge the gap between the two theories |^. For 
chaotic systems, each stable periodic orbit is locally near-integrable, and in general, will 
always be accompanied by island chains of non-isolated periodic orbits, which again inhibit 
the application of Gutzwiller's theory. To explore the limits of periodic orbit quantization for 
the case of stable motion we need to first investigate 'toy systems' of low generality like ours 
where the method is expected to work best. The circular resonances associated with each 
stable periodic orbit of our mapping rotate about each point of the orbit in a linear fashion, 
and hence, contain no high-order resonances. For almost all parameter values, the period of 
rotation will be an irrational multiple of tt. Thus, each stable periodic orbit will be isolated 
by the radius of its circular resonance, which, for a careful choice of the parameter values, 
can be made large enough such that Gutzwiller's theory is applicable. The low generality 
of our model should not render our research as being futile or purely academic. Indeed, the 
investigation of other toy systems, such as Sinai's billiard, has led to much insight for the 



hyperbolic case |31 



Our paper is sectioned as follows. In Section we review the classical map and construct 
its generating function. For a more detailed exposition of the classical map one should refer 
to |30[ . In Section |T| we introduce the quantum map and spin coherent states. Using these 



states as a basis, we then derive semiclassical matrix elements of the Floquet operator in 
Section |V|, followed by semiclassical traces and eigenphases in Section |V|. The theory in these 
sections follows the work by Kus et al. ||2^. Numerical investigations of our semiclassical 
approximations are contained in sections |VI| and |V11| . Finally, in Section Vlll we discuss 
our results. 
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II. CLASSICAL MAP 



The Hamiltonian under consideration is that of a kicked hnear top 

oo 

H{3,t) = ujj-i + n\Ji\ S{t-n), 



(1) 



where fi,uj E [0, 27r) are parameters and (J)j = Ji = eijkXjPk {i = 1,2,3) are the three 
components of angular momentum for a particle confined to a sphere, normalized such that 
J ■ J = 1. The classical evolution of J is governed by the equations 



Ji 



-ijkJk 



{Ji,H}, {Ji,Jj} 

where {■ , ■} are the Poisson brackets. Their solution can be written as the mapping 
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— sin yus" 
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(2) 



where s" = sgn J", which takes J from just before a kick to one period later. Here sgnx is 
the signum function with the convention sgnO = 0. 

The mapping is deceptively simple. It first rotates the eastern hemisphere (Ji > 0) 
through an angle /i and the western hemisphere (Ji < 0) through an angle — /i. The great 
circle Ji = remains fixed. Then the entire sphere is rotated about the J3 axis through 
an angle u. Thus the map rotates every point on the sphere in a piecewise linear fashion 
except on Ji = where its Jacobian is singular. Surprisingly, this leads to a phase space 
of great structural complexity. In Figure |l] we give examples for two different parameter 
values. The eastern hemisphere is shown in black, the western in gray. We will only consider 
the semiclassical quantization when fi = to = n{^/5 — 1) as in Figure |l](b). 

Note that if an orbit of our map does not have a point on Ji = 0, then its Lyapunov 
exponents are zero, making it stable. If however, it does have a point on Ji = 0, then 
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its Lyapunov exponents are undefined. We overcome this problem by simply defining an 
unstable orbit to be one with a point on Ji = 0, and all others stable. The unstable set [41] 
is defined to be the closure of the set of all images and preimages of the great circle Ji = 0. 
In figures |1] and |^ the unstable set is in black and gray. The circular holes are resonances 
consisting entirely of stable orbits. At their center lies a stable periodic orbit. It is these 
stable periodic orbits which will be used for the semiclassical quantization. All other points 
inside the resonances rotate about the central periodic orbit in a linear fashion. The period 
of rotation is an irrational multiple of tt for almost all values of the parameters /i and a;, 
and hence, no other periodic orbits are located inside the resonances. In our case of study 
(Figure |l|(b)), we have chosen /i to be large enough such that no island chains of first-order 
resonances form to create periodic quasi-tori (as in Figure |l|(a)). Hence we have done our 
best to ensure that all of the stable periodic orbits are sufficiently isolated. The unstable 
periodic orbits are shown not to affect leading terms in the asymptotic analysis. 



We have previously shown ||3^ that it is possible to label every point of a stable [42] 
periodic orbit of least period n with a unique sequence {s^ = ±1)^=1. The position of the 
point is given by the solution of 

J = F(s")F(s"-^)...F(s^)J = RJ (3) 

lying on the unit sphere with sgn Ji = s^. The other n — 1 points of the periodic orbit are 
found by substituting the n — 1 cycles of the original sequence into (|^). Hence, if {3^}k=i..n 
is a stable periodic orbit then the sequence corresponding to its first point is {sgn }fc=i..n- 
Note that the (unnormalized) solution of (^ is simply the axis of rotation of R, namely 

J = (R23 — R32 , R31 — Ri3 5 R12 — R2l)- 

Although every periodic point is uniquely represented by a sequence, not every sequence 
represents a periodic point. Hence we still need to find which sequences are legitimate. This 
may be done by checking each of the 2^ different possible sequences for a period-n orbit. 
However, this is computationally expensive for large periods. An alternative is to exploit 
symmetries of our mapping. We have previously conjectured [^ that every stable periodic 
orbit has exactly one or two points on one of the following symmetry lines on the unit sphere: 

Ji = —J2 tanco'/2 

J2 = J3tan(/i/2sgn Ji) 

J2 = Ji tanco'/2. 

It is easy to prove that if a stable periodic orbit has one point on a particular symmetry line, 
then it has no more than two points on this line. However, we were not able to show that a 
stable periodic orbit must have a point on one of these lines, though numerical investigations 
suggest they do. We have checked all 2^^ — 2 of the possible sequences for a stable periodic 
orbit of period < 40 and only found those of the above type. If one assumes our conjecture 
to be true then it is computationally easy to find all of the stable periodic orbits with periods 
upwards of 300 000. This may be done by finding where each symmetry line intersects with 
its image; an effortless task to solve numerically by virtue of the simplicity of our mapping. 
Every image of a symmetry line is just a collection of arc segments on the sphere. Once a 
periodic point is found, the size of its circular resonance can be calculated by noting that the 
resonance of at least one point of the orbit must touch Ji = 0. This is because the boundary 
of each resonance forms part of the unstable set, and hence, iterates arbitrarily close to 
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Ji = 0. The sections of the symmetry hne which intersect this periodic orbit's resonances 
contain no other periodic orbits and may be ignored for the remainder of our search. 

Following Kus et al. we rewrite the map in terms of the complex variables 7 and 7* 



7n+l 



In 



is'^ tan I 



i7„s" tan | 



r(7n)7n) 



(4) 



where 



7 



Ji + iJ2 



1 + 

is the stereographic projection of the unit sphere J 



e^'^tan 



fsin 9 cos d), sin 9 sin 6, cos 9) onto the 



complex plane. We then construct its generating function ( Re 7 7^ 0) 

Vr(r,7) =log [(l + 7)(l + e'"r) +e'^^(l-7)(l-e*"r*)l +C 



(5) 



where s = sgnRe7 [43] and C is an arbitrary constant. The mapping r(7,7*) can be 
rederived via the relations 



dW 



dW 



7 



ar* 1 + rr* 



d'j 1 + 77* 



(6) 



It has been shown that for nearly any area preserving map r(7, 7*) it is possible to find 
a generating function W{r*,'y) with 7 and T* as its independent variables with the above 
relations (^. Except at the discontinuity Re7 = 0, our particular map (^) is locally only a 
function of 7, F = F(7). Hence we also have the relations 



-7 



*2 



:i + FF*)2 



^7^ (1 + 77*) 



(7) 



III. QUANTUM MAP 

The solution of Schrodinger's equation with the Hamiltonian can be written in terms 
of a quantum map 

where F is a Floquet operator and the angular momentum operators Ji satisfy 

[Ji,Jj]=eijkJk, J^=j{j + 1) j = 0, |, 1, |, • • • 

For simplicity we have replaced J with hJ and noted that the normalization condition is 
now h^J'^ = 1 i.e. h = l/^Jj{j + 1). The classical limit is now approached by letting j 00. 
The quantum state lip"^) is a member of a Hilbert space spanned by the 2j + 1 orthonormal 
eigenstates of J3 

Mj^^) = m\j,m) 
( Ji ± z J2) I j, m) = ^J{j ^m){j ±m + l)\j, mil) 
where —j<m<j. By rewriting the Floquet operator as 
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where 

j 

m=-j 

one can find the matrix elements {j,m\F\j,n) and hence, the quantum traces trF"'. 

The nonorthogonal overcomplete set of spin coherent states make a more suitable 
basis for the semiclassical analysis. They are parametrized by the stereographic projection 
variable 7 and can be defined as a rotation of the minimum uncertainty eigenstate \j,j) [44] 

I7) = g*^('^lsin<^-J2Cos</>)|j 

= (1 + 77*)-^' E ir—^ ^V7^'+"^|J,"^) 

so that 

(7I JI7) = j(sin6' COS0, sin 9 sincj), cos9). 
The spin coherent states have minimal dispersion 

(7|AP|7) = (AP)^i„ = j 

and are thus closest to classical states. The trace of an operator in this basis 

tM = Hi±i/^!i,^(,iiw (9) 

IT J (1 + 77*)^ 



results from the following resolution of unity 

2j + l f cP-i 



TT 7(1 + 77*)' 

where (i^7 = (i( Re7)(i( Im7). 



■|7)(7l = l (10) 



IV. SEMICLASSICAL MATRIX ELEMENTS 



The semiclassical matrix element (r|F|7)sc is the leading term in the asymptotic expan- 
sion (j — i> 00) of the exact quantum matrix element. To proceed with its derivation, first 
consider the matrix element 



m=-j 

(2j)! 



j 

+ E 



(j + m) ! (j — m) ! 
m=o (j + m)!(j -m) 



(2j)! 



1+ 2Fi(-j,l;j + l;-eV- 



-in 



7 



where |^) and |?7) are arbitrary coherent states, and 2Fi{—j, l;j + 1; z) is a hypergeometric 
function. In Appendix ^ we derive an asymptotic expansion of this function when j is large. 
Using these results one obtains 



1 + 



(l + ee)^(l + r/r/ 



2j 



+ 



-4^*776 



1 + Cve''' 



if 1^*7/ 1 7^ 1, where s = sgn (1 — I^C*^!)- This approximation is exponentially small for large 
j except when ^ = e*^^//, in which case s = sgn(l — |?7|) and hence ^ is the image of r) 
under the classical map corresponding to e"*^'''^'. When this happens the first term in our 
expansion is 0(1) and exponentially dominates. Hence we take 

as the semiclassical approximation of the quantum matrix element. In the case |^*?7 
one finds that the leading term in the asymptotic expansion is 0(1) only when ^ = 



(12) 

I = 1 



with 1^1 = |?7| = 1. However the classical map ^ = e^'^^rj does not behave in this manner on 
the singular line |?7| = 1. In the next section we find that it is only when the mapping 1] ^ 
forms, through (^, part of a classical periodic orbit that contributions to the semiclassical 
trace become important. Hence we may assume |^*?7| 7^ 1, which amounts to ignoring the 



7-1 



unstable set of our 


original classical map 






Now using (|I2]) 


and the relations 






e'^^^|7) = 


(f 


+7y 








e-^-7) 



7 



with d^) one finds that 
where 



[1 + 7)(1 + e^'^r) + e^f^'il - 7)(1 - e^^F*) 



2i 



(13) 



s = sgn 1 



T - 1 



T + 1 



7-1 



In particular, if F = r(7, 7*) under the classical map 
matrix element can be rewritten in the two forms 



7 + 1 
then s 



(14) 

sgn Re 7. The semiclassical 



(r|F|7). 



dr*d-f 



exp 



(2j + i)w{T*, 7) (1 + rr)-^-(i + 77*)-^- 



exp 



2jw{r,j) (i + rr)-^a + 77*)- 



(15) 
(16) 



by putting C = — log 2 — iuj/2 — in/2 in (^. The first (p!5|) is reminiscent of that for the 
kicked top while the second (0) results from 



exp 



2W{r,j) 



(17) 



dT*dj 

which is special to this mapping. In both cases the semiclassical matrix element is fully 
determined by the classical generating function. 



8 



V. SEMICLASSICAL TRACES AND EIGENPHASES 



We now wish to derive the semiclassical approximation to the trace 



trF" 



2j + l 

TT 

2j + l 

TT 

2j + l 

TT 



:i+7i7r 



*\2 



(7il^l7i) 



(18) 



fc=i 

n 

fe=i 



+ 7fc7fc)' 
+ 7fc7fc)= 



:(7fe+il^l7fc) 



in + l 



exp 



2j(w^(7:+i,7.)-log(l + 7.7D)] (19) 



where we have used @, (|T0D and the semiclassical matrix element (|T6D. In the asymptotic 
limit j — > cxD this multiple integral is dominated by the contributions at saddle points where 



^W^(7fc+i,7fc) 



Ik+i 



k+l 



1 + 7fc+i7fc+i 



dWiYk+i,lk) 



1 + 7fc7fc 



for k = l...n{n+l = l). Recalling relations we see that the periodic points of the 
classical map are saddle points. However these are not the only saddle points. Contributions 
from 'ghost orbits' |^ also affect the above integral. Whenever a periodic orbit is destroyed 
at a bifurcation point there is still a residual contribution to the integral. These contributions 
are not semiclassical, being exponentially smaller, and will be ignored at this stage. However, 
near the bifurcation point the semiclassical approximation will become inaccurate if j is not 
large enough. By expanding the argument of the exponential up to second order one may 
approximate the contribution from each of the classical periodic points 



trF" 



E 



2j + l 



TT 



k=l 



■ exp 



2j exp 



2W^(Pfc+i,Pfc 



1 



7fe+i7fc 



1 



PkPk 



-Iklk 



(20) 



where p = {pk}k=i..n is the periodic orbit corresponding to the periodic point pi, and we 
have replaced 7 by 7 + p and used (0) and (|17D for simplification. This multiple integral may 
now be solved analytically (see Appendix ^). The result is our semiclassical approximation 
for the trace 



exp 


^(2j 


+ i)5(p); 


2i sin 


'S{p) 





(21) 



where the sum is taken over all periodic points (all unique points of each periodic orbit with 
period dividing n) and 

ts{p) ^ J2 wipi^,, p,) - iog(i + p^pi). 

k=l 



S is the action of the classical orbit 

Consider the contribution to the semiclassical trace formula from a single periodic orbit 
p = {pk}k=i..p with least period p and action S{p) from a single traversal. This orbit will 
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contribute when n = kp 

exp , 

(ti F ^)sc = p f 1 + (contributions from other orbits) 



exp 


i(2j + l)kSip) 


2i sin 


kSip) 





p ^ exp i{j — m)2kS{p) + (contributions from other orbits). 



m=0 

However the exact quantum mechanical trace is 

Ij-pkp = '^exp{-ikp(j)i) 

1=0 

where (pi are the eigenphases of the Floquet operator. Upon comparing these last two 
equations we see that for each m > there must be p distinct eigenphases with the property 

p(j)i 2{m — j)S{p) mod 27r as j — > oo, 

and hence, we obtain the following semiclassical approximation to the eigenphases 

(0m„p)sc = ^ -S{p) + mod 27r n = l,..p m = 0, 1, . . . (22) 

p p 

This approximation is not surprising. In the classical phase space each stable periodic point 
is at the center of a circular resonance, and all orbits in this resonance rotate about the 
periodic orbit in a linear fashion. Hence the map behaves, to a first-order approximation, 
like the linear top exp[2iS'(p) J3] near the periodic point. The point J3 = —1 of the linear 
top corresponds to the periodic point. When m = 0, the p different eigenphases correspond 
to the p different ground states associated with the periodic orbit. These ground states 
are localized on the periodic orbit and are p-fold 'cat states' (see Eq. (|27| ) and Fig. 
|^(a)). The higher order states are ranked via m and form rings around the classical periodic 
orbit similar to that of the linear top (see Figs. |^(b) and (c)). However, if j is finite, 
the semiclassical eigenphases will only be good approximates to the quantum for m ^ M, 
for some M, since only finitely many will 'fit' in the classical circular resonance. We can 
estimate M by considering the eigenstates of the linear top. These eigenstates have circular 
Husimi functions |(7|j, OP localized at an angle 6 satisfying 

tan^(^/2) = ^ -j<l<j. 

If we define Or to be the angular separation between a point of the periodic orbit and the 
boundary of its circular resonance, and find the eigenstate Ir of the linear top which is 
localized at this angle, then we obtain an estimate for M 



M{p) =j + lr = 2j sin^(^,,/2) = J - jy 1 - i?(p)2 (23) 

where 

R{p) = min \J^\ (= sin6',.)- 

k=l..p 

The calculation of R{p) follows from the fact that at least one of the circular resonances 
associated with the periodic orbit will be tangent to the great circle Ji = 0. 
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VI. NUMERICAL INVESTIGATIONS OF THE SEMICLASSICAL TRACES 



The semiclassical approximation to the traces pl| ) is exact in the hmit j — > cxd. But for 
finite j we need to test how large j must be before this approximation is an accurate one. 
We will only consider the case when fx = to = 7r{^/5 — 1). The actions S of all stable periodic 
orbits with period p < 25 have been compiled in Table |. Figure ^ shows the relative error 
of various semiclassical traces versus j 



,^trF"),c - trF^ 

relative error 

V|(trF«) 

where 



|2 



(24) 



;trF"),,|2 = ^ 



4sin^ 



S{P) 



Note that some of the traces can be quite slow to converge. In the case of n = 11 convergence 
is extremely slow. This was found to be caused by the failure of the semiclassical trace 
formula (0) when sin S'(p) =0. If n = 11 the only contributions to the trace are from 
two period-1 orbits each with action S mod 27r = 1.4391 from a single traversal. Hence 
sin 115* = —0.1218, which is close to zero making our Gaussian approximation (^Of ) an 
inaccurate one. Higher order terms from a complete asymptotic expansion of the integral 



(p!9|) would need to be calculated for increased accuracy. The problem is further exposed 
when the quantum diagonal matrix element (7|F"|7) is plotted. This is done for j = 300 in 
Figure H(a) using the stereographic projection 

Figure ||(c) shows the corresponding classical phase space. The periodic points of period 1, 
2 and 4 are labeled and the circle Ji = is shown in gray. One can see that an overwhelming 
contribution to the trace (p!8[) in the form of a large hump is produced by the presence of 
a classical period-1 orbit. However at j = 300 the hump does not take the form of a sharp 
Gaussian peak and is instead quite wide. 

Complete failure of the semiclassical trace formula occurs when sinS'(p) = 0. If this 
happens in the classical map the local rotation about the stable periodic point p, given by 



TABLE I: Act 


ions of all stable 


periodic orbits 


with period p <25 when /i = 


CO = 7r(v^-l). 


P 


S 


P 


S 


P 


S 


1 


-4.8441 


14 


-35.5088 


20 


-26.4359 


1 


1.4391 


14 


-10.3761 


20 


-39.0023 


2 


-0.6890 


16 


-20.0022 


20 


-64.5127 


4 


-6.4608 


16 


-45.1734 


20 


-1.6809 


9 


-17.7025 


16 


-7.4743 


22 


-48.6527 


9 


-11.2193 


18 


-54.8424 


22 


-23.5199 


12 


-13.2882 


18 


-4.5769 


24 


-32.8587 
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FIG. 2: The relative error (|2^) of various semiclassical traces versus j for jj, = uj = 7r(-\/5 — 1) 



2S{p), becomes a multiple of 2tt. In our case this occurs when fi = fii ^ u + 0.06556 where 
III is defined via 

1 + 2 cos j^n = cosuj + cos cos /ii + cos fii. 

The two period-1 orbits now have the action S mod 27r = 57r/ll. In more generic systems 
a period- 11 orbit would bifurcate from each period-1 orbit at this point. However no such 
bifurcation occurs in our mapping. Instead the circular resonance enclosing each of the 
period-1 orbits takes the form of an 11 sided polygon (Figure 0(d)). Every point inside this 
polygon rotates lOvr/ll radians about the period-1 orbit at each iteration of the mapping, 
and hence, is a period-11 orbit. The contribution to the semiclassical trace in these special 
cases is 

Mp) 



:trF"),, = (2j + l)5: 



where A{p) is the area on the unit sphere covered by the polygon. In our case we have a 
regular polygon with area 



A = 27T — 2n arcsin 



in ^y^l — Ji^ sin(7r/n)^ 



where n = 11 and the period-1 orbit is located at (Ji, J2, J3) with 



Ji = ±cos(u;/2)/^l + sin2(c^/2)cot2(/xi/2). 

In Figure ^ we have plotted the semiclassical and quantum traces for this case. The semi- 
classical displays a good first-order approximation to the quantum. Note that in the limit 
j ^ 00 the trace is infinite, whereas for = /^i ± e (e > 0) the trace limits to a large but 
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m 



FIG. 3: The quantum diagonal matrix element (7|F"|7) when j = 300 using the stereo- 
graphic projection (pSD, and with (a) n = 11, fi = u, (h) n = 11, fi = fii ^ u + 0.06556, (e) 
n = Q, n = oj, and (f) n = 6, /i = /ii. The stereographic projection of the classical phase 
space when (c) n = u and (d) /i = /ii. 
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FIG. 4: The semiclassical (black) and quantum (gray) trace trF^^ when /i = /zi. 



finite number. The semiclassical trace is this limit in the latter case. Consequently, when e 
becomes small (but still nonzero), the semiclassical trace ([21| ) will be an extremely poor ap- 
proximate for the quantum, since j needs to be ever larger before the quantum trace reaches 
its semiclassical limit. Hence the poor convergence in Figure when n = 11 is explained. 
Another test of the semiclassical traces is to consider the quantum Fourier transforms 

1 jo+N 

Tr,{9-jo,N) = - ^ e-'^HiF- (26) 

j=jo 

which exhibit peaks at ^ = 2S{p) mod 2tt. In Figures ^(a) and (c) we have plotted the 
absolute value of the quantum transform Tg when jo = 1 and jo = 1001, respectively, and 

= 200. Figure ^(e) shows the transform of the semiclassical trace. The peaks have 
been labeled according to the classical periodic orbit from which they derive. Note that 
the quantum transforms exhibit a third peak with no semiclassical analogue. This peak 
illustrates the effect of 'ghost orbits' [Q. It is the result of a residual contribution to the 
trace, left behind after a pair of period-6 orbits were destroyed when /j, uj + 0.05318. 
In Figure ^(f) we have plotted the quantum diagonal matrix element (7|F^|7) when fi = 
/ii ~ u; + 0.06556 and j = 300. One finds sharp Gaussian peaks at the location of the 
classical periodic orbits of periods 1, 2 and 6. The pair of period-6 orbits are marked by 
+'s and x's in the classical phase space (Figure 0(d)). In spite of there being no classical 
period-6 orbits when /i = cj, the quantum diagonal matrix element still displays ghost peaks 
at their previous locations (Figure 0(e)). These peaks, however, will vanish in the limit 
j oo unlike the case of /i = /zi where the peaks will instead increase to a height of unity. 
The exponential decrease in these ghost contributions is also exhibited in the quantum 
transforms (Figures ^( a) and (c)). Figures |^(b) and (d) show the quantum transform T25. 
The only contribution to the semiclassical trace is from the pair of period- 1 orbits (Figure 
0(f)). But in the quantum transforms this has been swamped by other contributions. The 
largest peak in Figure 0(d) is the ghost of a pair of period-25 orbits which were destroyed at 
fi ^ u + 0.0004288. This peak has temporarily grown in size after increasing j. The other 
quantum peaks are also believed to be caused by residual quantum effects. However we are 
currently unable to give an exact reason for their presence. 

The quantum transform T25 paints a disheartening picture of exponentially proliferating 
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FIG. 5: The quantum Fourier transforms Tn{9; jo, 200) when (a) n = 6, jo = 1, (b) n = 25, 
jo = 1, (c) n = 6, jo = 1001, and (d) n = 25, jo = 1001. The semiclassical transforms when 
(e) n = 6 and (f) n — 25. 



quantum effects not accounted for in the simple semiclassical analysis. In general we found 

our semiclassical trace (trF"')sc to be plagued by errors when n ^ 10. Although our 
approximation becomes exact in the limit j oo, and for large j the error is exponentially 
decreasing, this decay did not occur as quickly as we had hoped. In the next section we 
find that these inaccuracies make it impossible to extract semiclassical eigenphases from the 
traces. 



VII. NUMERICAL INVESTIGATIONS OF THE SEMICLASSICAL SPECTRUM 



Consider the spectral density of the eigenphases 
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FIG. 6: The semiclassical (black) and quantum (gray) spectral density p{9; j, N) when 
= 10 and (a) j = 50, (c) j = 100, (e) j = 1500, and when N = 30, (b) j = 50, (d) 
J = 100, (f) J = 1500. 



2j + l 

2n 



-J2 Re(e'"nrF"). 



TT 



n=l 



Hence, by rewriting the spectral density in terms of traces of the Floquet operator and using 
the semiclassical trace formula (pT]), one may derive a semiclassical spectral density. How- 
ever, this formula requires traces of the Floquet operator to all powers n, which invariably 
become inaccurate when n is large. Thus, we instead consider the spectral density with 
limited resolution 



piO;j,N) 



2j + l , 1 



N 



2tt 



Re(e™^trF") 



n=l 



In figures ||(a), (c) and (e) we plot the semiclassical and quantum (gray) spectral densities 
when = 10 and j = 50, 100 and 1500, respectively. Figures |^(b), (d) and (f) are the same 
except with = 30. The semiclassical approximation is of course most accurate when 
j = 1500 and = 10. However one needs A^ ^ 10 to resolve the 3001 eigenphases and even 
increasing to A^ = 30 one finds large discrepancies between the quantum and semiclassical 
spectral densities. 
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TABLE II: 


Semiclassical and quantum eij 


2;enphases 


when LL = uj = 


n(V5-l) and j = 100. 


m 


n 


P 


^ (Pmnp JSC 


W)qu 


A 
Z\ 





1 


1 


1.2064 


1.2064 


3 X 10^'' 





1 


1' 


1.2064 


1.2064 


4 X 10^'^ 


1 


1 


1 


4.0846 


4.0846 


2 X 10^® 


1 


1 


1' 


4.0846 


4.0846 


2 X 10^^ 


2 


1 


1 


0.6797 


0.6799 


2 X 10~^ 


2 


1 


1' 


0.6797 


0.6800 


3 X 10~^ 


6 


1 


1 


5.9093 


5.8913 


2 X 10~^ 


6 


1 


1' 


5.9093 


5.8533 


6 X 10~^ 





1 


2 


2.9306 


2.9306 


2 X 10 





2 


2 


6.0722 


6.0722 


1 X 10~^^ 


1 


1 


2 


2.2416 


2.2416 


1 in — 12 
1 X 10 


1 


2 


2 


5.3832 


5.3832 


1 X 10^^^ 


2 


1 


2 


1.5525 


1.5525 


1 X 10^^^ 


2 


2 


2 


4.6941 


4.6941 


1 X 10^^^ 


57 


1 


2 


1.3544 


1.3545 


1 X 10 


57 


2 


2 


4.4960 


4.4958 


2 X 10~^ 


58 


1 
± 


9 


0.6654 


0.6657 


3 X 10~^ 


58 


2 


2 


3.8070 


3.8006 


6 X 10-3 


59 


1 


2 


6.2595 


6.2590 


5 X 10^^ 


59 


2 


2 


3.1179 


3.1340 


2 X 10-2 





2 


4 


5.7403 


5.7420 


2 X 10-3 





4 


4 


2.5987 


2.6007 


2 X 10-3 


1 


4 


4 


5.6515 


5.6454 


6 X 10-3 



Another approach is to use Newton's formulae plf and rewrite the coefficients of 
the characteristic polynomial det(F — z) in terms of traces of the Floquet operator. All 
2j + 1 eigenphases (pm can then be extracted via the roots Zm = e-*'^" of the characteristic 
polynomial. However we still need accurate semiclassical approximates of the trace trF"' 
for 1 < n < 2j + 1 (or 1 < n < j + 1 if we enforce 'selfinversiveness' on the polynomial i.e. 
det(F — z) = z'^^'^^ det Fdet(F'^ — z^^) where det F = e~^'^^^^~^^^). The method of harmonic 
inversion 36 1 via filter-diagonalization |37, 38| also proved futile. The problem is not due 
to the extraction procedure but rather the inaccuracy of our semiclassical trace (trF")sc 
when n becomes large. 

Our semiclassical eigenphases (p2|), however, were found to provide good estimates of the 
quantum eigenphases. Some of the eigenphases for j = 100 are tabulated in Table The 
quantum eigenphases were calculated in the usual manner by diagonalizing the Floquet 
operator. The Husimi function |(7|^/')P of each eigenstate was then plotted in phase space 
to locate the periodic orbit from which it derives. Approximately 142 of the 201 quantum 
eigenphases could be classified as corresponding to a particular periodic orbit (120 for the 
period-2 orbit, 14 for the period-1 pair, and 8 for the period-4 orbit). 

In Figure |^ we have plotted the Husimi functions of selected eigenstates (refer to Figure 
H(c) for the classical phase space). One of the four different ground states associated with 
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m m 

FIG. 7: The Husimi functions of selected eigenstates when n — oj — 7r(-\/5 — 1) and j — 100. 
(a) (m,n,p)=(0,2,4), (b) (2,1,1), (c) (57,1,2), (d) (6,1,1'), (e) and (f) no association to a 
periodic orbit. 
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FIG. 8: The logarithm of the error in the semi classical eigenphases versus m for p = 1 (+) 
and p = 2 ( x ) . 



the period-4 orbit is plotted in Figure 0(a). The Husimi functions of the three other ground 
states are very similar to this one. In the limit j — oo, each of these states will be in the 
form of a 4-fold 'cat state' 



l^il7T/2\^^\ _|_ l^imn/2\ \ _|_ lgi7i7r/2 



|P4) 



(27) 



where k, I, m and n are integers, and \pk) is a coherent state centered at the /c-th iterate 
of the classical period-4 orbit. Hence, linear combinations of the four ground states with 
appropriate prefactors of the form e*^'^/^/2, will produce localized states on only one of the 
four points of the classical orbit. 

In Figure |^(b) we have plotted an eigenstate associated with the pair of period-1 orbits. 
This state is of higher order (m = 2) and forms rings around the two classical fixed points. 
The period-1 orbits have degenerate actions (mod 27r), and hence, their semiclassical eigen- 
phases are equal. We have primed the second of each quantum number p = 1 in Table |I| 
to emphasize this degeneracy. Note that the quantum eigenphases, however, are unequal. 



A small splitting allows 'quantum tunnelling' |39| between the period-1 orbits. If we denote 
the two ground eigenstates by |011) and |011'), then the state |011) + |011') will be localized 
on one of the period-1 orbits, while the state |011) — |011') will be localized on the other. 
Only after Tr/K^oii)^^ — (0oii')g«l = 3.05 x 10^ iterations of the quantum mapping, will each 
of these states have completely tunnelled their way across to each others period-1 point. 

In Figure |^(c) we have plotted an eigenstate associated with the period-2 orbit. This is 
another high-order eigenstate {m = 57), and despite the large value of m, it is immediately 
recognizable as being derived from the period-2 orbit. An example of when our association 
becomes questionable is given in Figure 0(d). We have classified this eigenstate as being 
derived from the period-1 pair. Figures 0(e) and (f) show examples of eigenstates which 
could not be classified as corresponding to a particular periodic orbit. Our semiclassical 
analysis affords no description of these eigenstates. 

The logarithm of the error in the semiclassical eigenphases 



A =1(0, 



mnp I sc 



qu I 



versus m is plotted in Figure ^ for the period-1 (+) and period-2 (x) orbits. This error 
decreases exponentially for decreasing m. The leveling out for m < 46 in the period-2 
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case is due to numerical error attributed to the use of double precision arithmetic. The 
dotted lines are our estimates for the maximum value of m (^) for which the semiclassical 
eigenphases will approximate a quantum eigenphase. In the period-1 case M = 5.98, whilst 
in the period-2 case M = 61.5. In both cases our estimates give reasonable approximations. 



VIII. CONCLUSION 

As a first attempt at applying Gutzwiller's method of periodic orbit quantization to 
the mapping under consideration we have had some success. However we have also met 
great difficulties. We have found that our semiclassical approximation for the trace tr F"' 
accumulates large errors upon increasing n. These errors decay for large j, however not at 
a rate rapid enough for resolution of the semiclassical spectrum. The errors were caused by 
an exponential proliferation of quantum effects, including ghosts, which were not accounted 
for in the simple semiclassical analysis. Despite this, our semiclassical eigenphase formula 
gave an accurate approximation to the quantum eigenphases associated with each periodic 
orbit. However, this accuracy is hardly surprising given the simplicity of our mapping in 
the neighborhoods of stable periodic orbits. The global dynamics is of more interest being 
of high complexity, similar to that of chaotic motion. The semiclassical analysis, however, 
could not describe eigenstates with support in these parts of the phase space. 

The accuracy of our semiclassical trace formula may be increased in three different ways: 

1. By including high-order terms in our approximation of the quantum matrix element 
([TT|) . These terms would need to be included if one were to study the semiclassical 
role of the unstable set. 



2. By including contributions from all saddle points of the integral (19). That is, include 
the exponentially small terms such as those from ghosts. 



3. By including high-order terms along with the first-order approximation pO| ) for the 
integral (|19]) . These terms would enhance the accuracy of our trace formula when the 
action S becomes close to a multiple of vr. 

However one must remember that, in general, an asymptotic series is expected to be diver- 
gent, and hence, including more terms may prove detrimental if j is not further increased. 

The usefulness of Gutzwiller's method of periodic orbit quantization is debatable. There is 
no doubt that the theoretical insights offered are second to none. However the practicality of 
applying such a theory to accurately calculate the eigenphases, or even the traces themselves, 
is extremely limited. One must always keep in mind that the semiclassical analysis is only 
valid in the limit j —>■ oo. How large j needs to be, before the theory becomes an accurate 
one, will depend on the specific nature of each particular system. In our case, Gutzwiller's 
method has proven inadequate. This failure may be related to the particular dynamical 
nature of our mapping, which is of extreme ellipticity rather than the preferred hyperbolicity. 
However we believe that the problems encountered in this study will, in general, reappear in 
other studies whenever large powers of the Floquet operator are considered. The theory is 
too simplistic to describe the global complexity encountered in most chaotic systems. Only 
when j is made large enough, such that some eigenstates will be dominated by the behavior 
in the immediate vicinity of a periodic orbit, will the theory produce accurate results. This 
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view may be a depressing one, however it is hardly surprising given the particular nature of 
the Gutzwiller's approach. 

In conclusion, we are left with yet another, albeit cute, example of periodic orbit quan- 
tization in an ever expanding sea of literature on the subject. However we hope that our 
contribution is judged favorably, being qualitatively different from the much explored hy- 
perbolic periodic orbit quantization. 
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APPENDIX A: ASYMPTOTICS OF 2Fi{-j,l] j + I; z) 

The hypergeometric function under consideration is first replaced by another using a 
recurrence relation [45] 

1 1 

2Fi(-j,l;j + l;2;) = - + ^ (l - j, 1; j + 1; ^) . (Al) 

This new hypergeometric function can be represented by the integral 

2Fi(l - J, l-j + l-z)=jf^[{l- t){l - zt)y-'dt. (A2) 

When |z| < 1 a complete asymptotic expansion can be found with the help of a quadratic 
transformation [46] 

,F,(l - j, l;j + l;z) = l;j + 1; ^jf^'j (A3) 

- 1 gl02±iZMf^)\o,r"). (A4) 



l + z^^ 0F(j+m)! {{l + zf 

Taking the limit n oo the series only converges when |42;/(1 + 2;)^| < 1, however we 
still believe it gives a valid asymptotic expansion of (|A2|) if \z\ < 1. This is supported by 
considering a path of steepest descent in the integral. Letting 

fXt) = {l~t){l-zt) 

we see that an asymptotic expansion of (|A2|) will involve a contribution from the endpoint 
t = where f(t) = 1, and perhaps a contribution from the saddle point t = 1/2(1 + 1/ z) 
where f'{t) = 0. Letting t = x+iy, the path of steepest descent from t = is given implicitly 
by Im (/(x + iy)) = 0. When \z\ < 1 this path is 



y{x) = 1^1+ Rez -2x - ^{1+ Rez - 2xf + 4(1 - \z\^)x(l - x) 

which connects t = to the other endpoint t =1 where /(t) = 0. Along the contour /(t) is 
monotonic, and hence, the expansion depends only on the local behaviour near t = 0. On 
deriving the first terms in the expansion one finds agreement with (|A4j) . 
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When > 1 the path of steepest descent from t = does not pass through t = 1. It 
does, however, pass through the other zero at t = To investigate this case we start by 
rewriting as 

2Fi(l-j,l;j + l;z) =j f'^"[(l-t){l-zt)y-'dt + j f [{I - t){l - zt)Y-^dt 

■/ J \ I z 

= ^- - s){i - s/z)y-^ds + ji-z)'^{i - zf^~^ [\s{i - s)y-^ds 

z Jo Jo 

where we have substituted t = s/z in the first integral and t = l/z + s(l — l/z) in the 
second (s real). The contour in the second integral takes a path of steepest ascent from 
t = 1/2; to the saddle point and then of steepest descent to t = 1. This integral can be 
evaluated explicitly, giving an exact contribution from the saddle point. The contour in the 
first integral is not of steepest descent but now this integral takes the form of (^). We have 
in fact rederived the linear transformation [47] 



2Fi{l-j,l;j + l;z 



2F,(l-j,l;j + l;l/z)+2y ^ 



z 



(2j)! 



'1-z 



(A5) 



The asymptotic expansion of ( |A^) when |2;| > 1 can now be found by using (|A4|) . 

In the case 1^1 = 1 {z 7^ —1), the path of steepest descent is a straight line from t = to 
the saddle point and then to t = 1. Hence we rewrite (^3) as 



- J, l;j + l;z)=j ' ' [(1 - t)(l - zt)r^dt^3 I J(l " - zt)r' 



'dt 



j(l + l/z) / ' [1 - 4cos2(e/2)s(l - s)Y~'ds + 



(jO! 

(2j)! 



:i-z 



,2i-l 



where z = e*^ and we have substituted t = s(l + l/2;)in the first integral and evaluated the 
second. When z ^ —1 the first integral takes its maximum at the endpoint s = 0. Again, the 
contribution from this point is (|A4|) . There is also a contribution from the endpoint s = 1/2. 
This contribution is exponentially smaller, being 0(j~") for all n, and will be ignored. If, 
however, j is not large enough the resulting expansion will become inaccurate when z is close 
to —1. Hence we take ( [A^ ) with the result of second integral as our asymptotic expansion 
in this case. 

Summarizing the above results, when z ^ —1 the asymptotic expansion of our original 
hypergeometric function is 



1 



1 — sgn (1 — \z\) 



'l-z) 



— z 



+ 



+ 



11 - ;z 



n— 1 

E 



r(m+ l/2)j! 



Az 



and when z 



21 + z^o \il + z)\ 

-1 we have the exact result from (|A5|) 



m 2- 



+ 0{rl (A6) 



(A7) 



We have tested these asymptotics numerically and found good agreement. However conver- 
gence is weak near z = —1. 



22 



APPENDIX B: A MULTIPLE INTEGRAL 

Consider the relation 



„ n 

In{a^,...,aJ = J Y[d 



k=i 1 — ai . . . 



where we have assumed the integral is convergent, z = x + iy, d?z = dxdy and the subscript 
n + 1 = 1. The proof is relatively simple and is achieved through induction. When n = 1 
the integral is trivial. We now assume the equality holds for n = m and consider the case 
n = m + 1 



J ,Ja. a -,1= d'^Z. d'^Z , p"'l4%+---+"m<+l2m+«m + l4^m+l-4^1---^m + l^m + l 

= 71 d'^ Z-, Cp Z g'^1^2^1+--- + '^m-l^m^m-l+'^m'^m+l^r^m~^iL^l^---^^m^m 

J 1 • • • m 



.m+1 



7T 



1 — ai . . . a^a. 



Hence the equality is true for all n. By making appropriate substitutions for and z^ the 
multiple integral (pOD is solved. 
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